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Abstract 

We perform an extensive study of the DAMA annual modulation data in the 
context of inelastic dark matter. We find that inelastic dark matter with mass 
m x > 15GeV is excluded at the 95% confidence level by the combination of DAMA 
spectral information and results from other direct detection experiments. However, 
at smaller m x , inelastic dark matter constitutes a possible solution to the DAMA 
puzzle. 
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1 Introduction 



The DAMA collaboration has recently published the combined results from the long- 
term runs of the DAMA/Nal and DAM A/LIBRA apparatus. Both experiments observe 
an annual modulation in their nuclear recoil rates with a combined statistical significance 
exceeding 8a [1]. A possible origin of this signal is weakly interacting massive particles 
(WIMPs) scattering off of the target nuclei in the detector, where the modulation of the 
rate is caused by the motion of the earth around the sun. The amplitude, period and 
phase of the signal are all consistent with such scatterings. 

This interpretation is challenged, however, by the absence of a signal in other dark 
matter direct detection experiments like CDMS [2] and XENON [3]. Especially in the 
standard scenario where WIMPs scatter elastically and spin-independently off of the 
target nucleus, the DAMA preferred region of WIMP cross sections is practically ruled 
out by the combination of other experiments [4,5]. 

In order to reconcile the DAMA signal with the other direct detection experiments, 
a number of non-standard dark matter candidates have been proposed. These include 
mirror dark matter [6,7], dipolar dark matter [8], hidden charged dark matter [9] and 
composite dark matter [10] as well as particles which couple spin-dependently [11] or 
preferentially to electrons [12,13]. 

In the present article we focus on inelastic dark matter [14] which is one of the 
simplest possible solutions to the DAMA puzzle. Concrete realizations of inelastic dark 
matter in particle physics can e.g. be found in [14-17]. In this scenario the WIMP scatters 
from a ground state x into the slightly heavier state x' ■ If the mass splitting 5 = m x i — m x 
is of the same order as the kinetic energy of WIMPs in the galactic halo, only the WIMPs 
in the tail of the velocity distribution are energetic enough to scatter off of a nucleus. In 
this case the the fraction of WIMPs which can overcome the kinematical barrier depends 
very strongly on the velocity of the earth, and the annual modulation of the WIMP 
scattering signal is significantly enhanced. This raises the sensitivity of DAMA relative 
to other experiments and reduces the tension between them. Earlier studies indeed found 
parameter regions consistent with DAMA and all other experiments for WIMP masses 
m x > 50 GeV and 5 = 0(100 keV) [16, 18, 19]. 

In this article we improve on those earlier studies. In particular we take into account 
the energy resolutions of the DAMA/LIBRA and DAMA/Nal apparatus and include 
higher energy bins which are typically ignored. We also use more appropriate analysis 
techniques than some of the earlier works. Furthermore we extend the analysis to small 
WIMP masses and find a region at m x ~ 0(10 GeV) where the DAMA signal can be 
explained through channeled scattering events off of iodine. Another allowed region at 
m x ~ 5 GeV - where the modulation is due to channeled sodium events - opens up if 
one artificially decreases the detector resolutions of DAMA/LIBRA and DAMA/Nal. 

We then confront the DAMA allowed regions with the null results of various other 
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direct detection experiments. In contrast to previous studies we find that under standard 
astrophysical assumptions the whole high mass region is ruled out at the 95% confidence 
level (CL) by CRESST II. For this statement to hold it is important that we do not 
combine the CRESST II commissioning run [20] with an older test run [21], as will be 
discussed in more detail later. 

The low mass regions corresponding to channeled scattering events at DAMA are 
partly probed by low-threshold experiments. However, here we still find regions at 
m x ~ 15 GeV which are compatible with DAMA as well as all other direct detection 
experiments. 

The outline of this paper is as follows: in the next section we review relevant aspects 
for the direct detection of inelastic dark matter. In Section [3] we give an overview over 
the direct detection experiments we consider in this study and describe our analysis 
procedure. In Section [4] we present the results of this analysis which we compare with 
results from previous studies in Sectional Finally Section [6] contains our conclusions. 

2 Direct Detection of Inelastic Dark Matter 

Direct detection experiments aim to measure the deposited energy of a WIMP dark 
matter particle when it interacts with a nucleus in the detector. If such a dark matter 
particle x scatters off of nuclei inelastically by making a transition to a slightly heavier 
state x'i the minimum incident velocity to transfer the energy Er to a recoiling nucleus 



where 5 is the mass splitting between x an d x', fn N is the mass of the target nucleus and 
fi = {m x m^) / {m x + m^) is the reduced mass of the WIMP-nucleus system. The altered 
kinematics of the scattering process relative to elastic scattering can affect the detection 
rates significantly. We will discuss these effects in more detail after a short reminder of 
the relevant formulae for the calculation of reaction rates and a brief discussion of the 
relevant input parameters. We will also specify the impact of detector effects on reaction 
rates. 

2.1 Reaction Rates 

The differential event rate for coherent spin-independent WIMP nucleus scattering as a 
function of the recoil energy E R can be written as (see e.g. [22]) 



is [14] 




(1) 



oo 



dE R 



dR 



(E R , t) = M x 



tar 




/local (V,t) 



V 




with M tar the mass of the target detector, m x the WIMP mass and p x the local WIMP 
density. Furthermore f n>p are the effective coherent couplings to the neutron and proton 
respectively, while A and Z denote the nucleon and proton numbers of the target nucleus 
and a n the overall effective WIMP-neutron cross section at zero momentum in the elastic 
limit. The nuclear form factor F(E R ) describes the loss of coherence as the momentum 
transfer deviates from zero. Finally, the integral takes into account the effect of the local 
WIMP velocity distribution /i oca i- Both, the input from particle physics as well as the 
astrophysical input together with the corresponding uncertainties will be addressed more 
closely below. 

2.2 Input Parameters 

To keep our analysis as general as possible, we keep the WIMP mass m x , the mass 
splitting S and the cross section a n as free parameters. These parameters may be predicted 
by the particle physics model under consideration. 

2.2.1 Input from Particle Physics 

For this study we set the ratio of the effective WIMP-proton and WIMP-neutron cou- 
plings to one, f p / f n = 1. This ratio could be different in a specific model. However, 
a ratio f p /f n ^ 1 can easily be corrected for. One just has to rescale the considered 
constraint on a n by a factor 



As all nuclei have roughly the same ratio of protons to neutrons, a modification of f p / f n 
does virtually not change the experimental constraints relative to each other. 

In our analysis we use the Fourier-Bessel form factors where available, otherwise 
the Woods-Saxon form factors^. Both can be calculated from the parameters tabulated 
in [23,24]. Note that Fourier-Bessel and Woods-Saxon form factors are more accurate 
than the commonly used Helm form factors (see e.g. [25]). Especially for heavy elements 
the Woods-Saxon form factors can deviate from the Helm form factors substantially. For 
example in the case of a tungsten target as used by CRESST II the Helm form factor 
would lead to an overestimation of the reaction rate by roughly 10%. 

2.2.2 Astrophysical Input 

The two main astrophysical inputs used in this study are the local dark matter density 
p x and the dark matter velocity distribution seen by an earth bound detector. For our 

1 Fourier-Bessel form factors are available for oxygen, aluminum, silicon and germanium, Woods- 
Saxon form factors are used for sodium, iodine, xenon and tungsten. 




(3) 
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analysis, we set p x = 0.3 GeVcm 3 . Note that a change of p x would just lead to an 
overall shift of all experimental constraints on a n . 

The galactic dark matter halo is expected to have a rather smooth velocity distri- 
bution. The Standard Halo Model, which we adopt throughout this study, assumes a 
Maxwell-Boltzmann distribution in the galactic rest frame, which is truncated at some 
escape velocity t> esc [26], 

= (^c - H) • (4) 

The local galactic escape velocity was determined to be in the interval 498 km/s < v esc < 
608km/s at 90% CL by the RAVE survey [27]. We will use the median v esc = 544 km/s. 
The velocity dispersion is set to its standard value v = 220 km/s. 

The dark matter velocity distribution in the rest frame of the earth, /i 0C ai 5 can be 
obtained from the galactic distribution by a Galilean boost, 

/local (V, t) = f gal (v + V earth (t)) , (5) 

where the velocity of the earth relative to the galactic halo can be parameterized as 



dearth (t) — V sun + faring) • 

Using the convention for galactic coordinates as in [28] we have 



km/s 



describing the motion of the solar system [29,30] and 
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corresponding to the time-dependent velocity of the earth relative to the sun [22,28]. Here 
t is measured in years and t\ = 0.219 corresponds to the Spring equinox on March 21. 
The maximal and minimal velocities of the earth with respect to the galactic rest frame 
are reached on June 2 and December 2 respectively. 



2.3 Annual Modulation and Inelastic Dark Matter 

The motion of the earth around the sun leads to an annual modulation of the WIMP 
velocity distribution with respect to the earth. Due to this annual modulation, the dif- 
ferential rate of WIMP scattering has a constant and a nearly sinusoidal time-dependent 
contribution, 

d-R 

—(E R ,t)=S + S m coB(2<ir(t-t )), (9) 
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with t corresponding to June 2nd, where the relative earth- WIMP velocity is maxi- 
mal. Note that for non-standard velocity distributions, which have additional features, 
the time dependence could in principle have a more complicated form. The modulation 
amplitude of the signal S m is given by 

S m (E R ) = \ June 2) - ^(E R , Dec 2)) . (10) 

Since DAMA, unlike other direct detection experiments, does not attempt to suppress 
backgrounds completely, it can only identify the modulated part of the signal. Taken 
by itself, the DAMA signal is consistent with the standard scenario of elastic WIMP 
scattering, even though the modulation S m / So is expected to be only at the level of a 
few percent in this case. This interpretation is, however, inconsistent with several other 
direct detection experiments and if the DAMA signal really is due to WIMP-nucleus 
scatterings, other dark matter interpretations have to be considered. 

One proposal to reconcile the DAMA signal with the null results from other ex- 
periments is inelastic dark matter. The inelasticity of the scattering drastically changes 
the kinematics of the process, c.f. (ED). In particular the higher kinematical barrier for 
inelastic processes allows only WIMPs in the tail of the velocity distribution to scatter. 
This leads to a strong enhancement of the annual modulation of the scattering rate - in 
the extreme case that the kinetic energies of WIMPs needed to scatter inelastically are 
reached only during summer, the modulation can even be maximal! This effect raises 
the sensitivity of DAMA compared to other direct detection experiments and ameliorates 
the tension between them. Furthermore, as the minimal velocity (DQ) typically decreases 
at higher nucleus mass, experiments with heavy target nuclei can detect a larger frac- 
tion of the WIMPs passing through the detector. This implies e.g. that the sensitivity 
of CDMS is reduced relative to DAMA, since both germanium and silicon are consid- 
erably lighter than iodine. On the other hand xenon is similar in mass and hence the 
XENON10 experiment should be sensitive to inelastic recoils. But since XENON10 took 
data only during winter where the reaction rates are smallest, the resulting limits can 
still be evaded. Nevertheless the other experiments still very strongly constrain the in- 
elastic dark matter interpretation, with especially strong constraints arising from the 
CRESST II experiment, which has tungsten as the heaviest target nucleus. 

2.4 Detector Effects 

In a real experiment there are additional effects which have to be taken into account to 
understand a possible signal. First of all, only a fraction of the energy of the recoiling 
nucleus will be transferred to the channel which is observed, typically as ionization or 
scintillation in the detector. This fraction is known as the quenching factor Q which 

2 See e.g. the left panel of Figure [2] for the typical annual modulation of the recoil rates. 
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relates the recoil energy E R to the observed energy E' through 

E' = QE R . (11) 

The observed energy E' is often referred to as electron equivalent energy and measured in 
keVee. The quenching factor mainly depends on the target material of the detector. Some 
experiments calibrate their energy scales such that the quenching factor is effectively one. 
For experiments with Q 7^ 1 one can apply (TTTT) to obtain the measured events in terms 
of the recoil energy. 

Assuming an ideal experimental setup the event rate would simply be given by (J2J). 
However, a realistic detector has only a finite energy resolution, i.e. for a recoiling nucleus 
with energy E R the measured energy E' follows a probability distribution which is peaked 
at Q Er. Typically one approximates the distribution by a Gaussian with an energy 
dependent standard deviation c(E'). To obtain the predicted differential reaction rate 
corrected by the finite energy resolution of the detector, one must convolute dR/dE R 
from (J2J) with this Gaussian, 

00 

dfl/r, x 1 / in/ dfl /n , , 1 / (E R -E') 2 \ . . 

jpr(E R ,t) = -= / dE' R -—(E' R ,t) — — exp - \* f * , 12 

^Corrected V^TT J dE R (j{E R ) \ 2a 2 (E R ) J 



where cr(E R ) may be gained from cr(E') by use of the quenching factor. 

Finally there is an energy dependent efficiency £, which takes into account the loss 
of signal, e.g. due to data cuts which are needed to reduce background. To determine the 
number of recoils which are expected to be seen at a given experiment between recoil 
energies E\ and E 2 in a time interval ti-t 2 , we have to integrate over the resolution 
corrected rate multiplied by the efficiency, 

ti E 2 

N El _ E2 = [dt [ dE R £(E R ) 4^-( E Rit) ■ (13) 



dE 



R Corrected 



For detectors with multiple elements/isotopes, the contributions for each element /isotope 
have to be added. 



3 Direct Detection Experiments 

In this section we discuss key properties and results from the DAMA experiment as 
well as from various null searches. We restrict ourselves to experiments which provide 
sufficient data to independently generate constraints. 

3.1 DAMA 

The DAMA collaboration has recently published the observed modulation amplitudes in 
36 energy bins over 2-20 keVee, where the data from DAM A/LIBRA and DAMA/Nal 
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were combined. We will shortly describe how to calculate the modulation amplitude 
expected at DAMA and then outline the statistical procedures used to determine the 
confidence regions for DAMA. 



3.1.1 Modulation Amplitude at DAMA 



DAMA's target consists of highly radiopure Nal(Tl) crystals. In principle both compo- 
nents - sodium and iodine - may contribute to the observed signal. Their respective 
quenching factors were measured to be Qn& = 0.3 and Q\ = 0.09 [31] respectively. 

While typically only this small fraction of the recoil energy goes into scintillation 
which can be measured by DAMA, some nuclei that recoil along the characteristic axes or 
planes of the crystal structure may transfer their energy primarily to electrons resulting 
in Q ~ 1 [32]. This effect is known as channeling. DAMA has used simulations to 
determine the energy-dependent fractions of channeled events off of sodium and iodine. 
These fractions can be extracted from Figure 4 in [32]. A good analytical approximation 
is given by 



0.63 R V 1+ i+^, 



and 



R 



J70.45 / 

0.66 R 1 



1+E R 



(14) 



respectively, where the recoil energy E R is understood in units of keV. In order to be able 
to compare predicted and measured rates, it is useful to calculate the predicted spectrum 
in terms of the measured energy. Taking into account quenching and channeling and 
adding contributions from sodium and iodine, one finds 
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(15) 



where dR\/dE R and dR^ & /dE R may be gained from j2]), while 0.847 and 0.153 denote 
the mass fractions of iodine and sodium in Nal respectively. The event rate has then to 
be corrected by the detector resolution according to ( fl~2l) . For DAM A/LIBRA cr(E') was 
determined to be [33] 



^dama/libra(£') = (0.448 keVee) 



E' 



keVec 



+ 0.0091E' 



(16) 



The resolution of DAMA/Nal is slightly worse, it can be obtained by a simple fit to 
Figure 13 in [34], 



o"DAMA/Nai(£') = (0.76 keVee) 



E' 



kcVcc 



0.024E'. 



(17) 
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The different energy resolutions of the two experiments force us to calculate the resolution 
corrected spectra of DAMA/LIBRA and DAMA/Nal separately. We then combine them 
weighted by their relative exposures which are 65% and 35% respectively, 



^-,0=0.66 &>,<)) + 0.36 

U-C/ Corrected \ U-C/ Corrected/ DAMA/LIBRA \ U-C/ Corrected/ DAMA/Nal 

Finally the predicted modulation amplitude averaged over an energy bin reads 



(18) 



*2 

*d = 2fF / F>] [ ^ (^(E', June 2) - ^(E>, Dec 2)) . (19) 

4 — ^l) J V Corrected Cl-C/ Corrected / 



3.1.2 Calculating Limits 

To find regions in the parameter space of inelastic dark matter which are consistent with 
DAMA, we employ a x 2 -goodness-of-fit metric, 

( cbin cbin \ 2 

2 _ V pred dataJ (r>p , 

bins \ u data) 

where 5"^° a stands for the measured modulation amplitude, while cr^ 1 ^ denotes the 
corresponding experimental error. The modulation signal is consistent with zero above 
10 keVee, therefore it is not appropriate to simply fit to all 36 energy bins of DAMA. 
Fitting to a large number of bins, where the measured modulation amplitude fluctuates 
around zero with a relatively large error, would just dilute the power of the goodness 
of fit metric. In [16, 18, 19] this problem was circumvented by just ignoring the higher 
energy bins of DAMA. However, we follow [5] where all bins above lOkeVee are combined 
into a single bin. The inclusion of this additional bin makes sure that scenarios in which 
the modulation amplitude is non-zero above 10 keVee are strongly constrained. 

In total we use 17 energy bins for the goodness-of-fit test, cf. Table [TJ As we are fitting 
to 17 energy bins, the 90% and 99%-confidence levels are characterized by x 2 < 24.8 and 
X 2 < 33.4 respectively. 

The DAMA collaboration has also published the absolute rate So (cf. (J9j)) measured 
by the DAMA/LIBRA apparatus in the range 0.75-10 keVee (cf. Figure 1 in [1]). As 
there is only limited background suppression, a large fraction of the measured rate may be 
due to background events. However, this data set can still be used to derive conservative 
exclusion limits bv imposing that the predicted absolute rate may at least not exceed 
the measured ratqj. For this constraint we use only the energy bins above 2 keVee as the 
behavior of the detector below its threshold may not be well understood. 



3 One does not have to include a statistical error in this calculation as the number of events observed 
in each bin is extremely large. 
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Energy 


Cbin 
°data 


Energy 


Cbin 
°data 


[keVee] 


[cpd/kg 


/keV] 


[keVee] 


[cpd/kg/keV] 


2.0- 


2.5 


0.0162 ± 


0.0048 


6.5-7.0 


0.0012 ± 0.0036 


2.5- 


3.0 


0.0287 ± 


0.0054 


7.0-7.5 


-0.0002 ± 0.0036 


3.0- 


3.5 


0.0250 ± 


0.0055 


7.5-8.0 


0.0004 ± 0.0036 


3.5- 


4.0 


0.0141 ± 


0.0050 


8.0-8.5 


-0.0014 ± 0.0037 


4.0- 


4.5 


0.0100 ± 


0.0045 


8.5-9.0 


0.0039 ± 0.0037 


4.5- 


5.0 


0.0118 ± 


0.0040 


9.0-9.5 


-0.0033 ± 0.0037 


5.0- 


5.5 


0.0039 ± 


0.0040 


9.5-10.0 


-0.0070 ± 0.0038 


5.5- 


6.0 


0.0030 ± 


0.0039 






6.0 


6.5 


0.0060 ± 


0.0038 


10.0-20.0 


0.0010 ± 0.0008 



Table 1: DAMA modulation amplitudes (and errors) in the lowest sixteen energy bins between 
2 and 10 keVee, extracted from Figure 9 in [1]. Bins above 10 keVee are combined into a single 
bin. 

3.2 Null Searches 

Here we list some of the experiments which give the most stringent bounds on the dark 
matter interpretation of the DAMA signal. Key parameters of the experiments we use 
can be found in Table El While most of these experiments observe some recoils, these 
recoils may be due to background events. In this work we do not attempt to subtract 
the background, hence we obtain conservative constraints. 

For unbinned data, constraints are generated using S. Yellin's maximum gap method 
[40], which examines the likelihood of observing the gaps in energy between observed 
events. This unbinned method generally provides a stronger constraint than binned 
methods for the case of an unknown background. Unfortunately experimental data is 
often available in binned form only. In this case we use the 'binned Poisson' technique, 
as described in [5,41]. 

3.2.1 XENON10 

The XENON10 detector acquired 316.4 kg-days of data between October 6, 2006 and 
February 14, 2007 and observed 10 candidate events [3]. Here energies are calibrated to 
recoil energies, implying Q = 1. This calibration depends upon the relative scintillation 
efficiency £ e fj which was assumed to be constant, £ e ff = 0.19 [3]. However, new measure- 
ments [42,43] indicate that the scintillation efficiency decreases to £ c g ~ 0.14 for recoil 
energies Er < 10 keV (see Figure 8 in [43]). We account for the lower C c r by rescaling 
the lower threshold of XENON10 and the lowest energy event by a factor 0.19/0.14. All 
other events are located at larger recoil energies and are not affected. After our rescaling 
the energy range of XENON10 is 6.1-26.9 keV. 



10 



Experiment 


XENON10 


ZEPLIN III 


CDMS II 


CDMS-SUF 


Target 


Xe 


Xe 


Ge 


Ge 


Exp. [kg-day] 


316.4 


126.7 


397.8 


65.8 


Range [keVee] 


6.1-26.9 


2-16 


10-100 


5-100 


Q 


1 


P5 


1 


1 


Efficiency [Eq.] 


m 




(25)) 




Resolution [Eq.] 


(E2D 




(26j) 


(26} 


h[y] 


0.763 


0.158 


~ 0.80 


0.995 


h[y] 


1.122 


0.382 


~ 1.55 


1.248 


Constraint 


maximum gap 


maximum gap 


maximum gap 


maximum gap 


Experiment 


CRESST I 


CRESST II 


CoGeNT 


TEXONO 


Target 


Al, 


W 


Ge 


Ge 


Exp. [kg-day] 


0.80, 0.71 


30.6 


8.38 


0.338 


Range [keVee] 


0.6-20 


10-100 


0.34-4.1 


0.2-7.1 


Q 


1 


1 


0.2 


0.2 


Efficiency [Eq.] 


1 


0.9 




ii 


Resolution [Eq.] 


(EHD 


(Ei 




(33} 


h[y] 


~ 0.75 


0.234 


~ 0.17 


0.040 


h[y] 


~ 0.75 


0.558 


~ 0.25 


0.100 


Constraint 


binned Poisson 


maximum gap 


binned Poisson 


binned Poisson 



Table 2: Parameters of direct dark matter searches for XENON10 [3], ZEPLIN III [35], CDMS- 
SUF [36], CDMS II [2], CRESST I [37], CRESST II [20], CoGeNT [38] and TEXONO [39]. Note 
that for CoGeNT and Texono we only use the energy bins below 2 keVee in our analysis (see 
text). 
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The energy dependent efficiency can be approximated by (c.f. Table 1 in [3]) 
£(E R ) = 0.46 f 1 - -Jfe- ) (21) 



100 keV, 

while the energy resolution can be written as [5] 



a(E R ) = (0.579 keV)^^| + 0.021£ fl . (22) 
3.2.2 ZEPLIN III 

Between February 27 and May 20, 2008, ZEPLIN 111 acquired 126.7 kg-days of data [35], 
taking into account all the energy independent efficiencies. Seven events were detected in 
the observed energy range 2-16 keVee, corresponding to recoil energies 10.7-30.2 keV, 
cf. Figure 12 in [35]. The quenching factor Q has to be extracted from the conversion 
between the energy scales in Figure 15, an analytic approximation is given by [19] 

Q(E') = (0.U2E' + 0.005) Exp(-0.305£' a564 ) . (23) 

The energy dependent efficiency is also shown in Figure 15 of [35]. In the considered 
window of recoil energies it is well approximated by 

. nfr1 15.23 keV 197.3 keV 2 , „ 

£(E R ) = 0.71 + — — 2 . (24) 

Unfortunately there is no information available on the energy resolution of ZEPLIN III. 
However, as XENON10 and ZEPLIN III have a similar detector setup, we will em- 
ploy (ED) also for ZEPLIN III. 



3.2.3 CDMS II 



CDMS II took data from October 2006 to July 2007, resulting in an exposure of 397.8 kg- 
days on its germanium detectors with an energy threshold of 10 keV [2]. No events were 
observed. The quenching factor is Q = 1. The efficiency of observing nuclear recoils can 
be extracted from Figure 2 in [2] (arXiv version only), we use the approximation 

f 0.24 + - 35 ' 3 Jf v2 for 10 keV < E R < 40 keV , 

Z(E R ) = l Er (25) 

I 0.30 for 40 keV <E R < 100 keV . 

The energy resolution for the CDMS germanium detectors is given by [44] 

a{E R ) ~ 0.2^ • (26) 
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3.2.4 CDMS-SUF 



Due to the low 5keV threshold we also include the 65.8 kg-days germanium data set from 
CDMS-SUF [36] in our analysis, which was taken from December 2001 to April 2002. 

The total efficiency can be approximated by [2] 

, s fo.lO + 0.30 ( ^~ 5 t eV) for 5 keV < E R < 20 keV, , x 

£(E R ) — 0.80 x 0.95 x < fp 15 ^ n 27 

1 0.40 + 0.10 ( fl 8 Q ke v ^ 20 keV < E R < 100 keV. 

In total there were 20 candidate events in the germanium detectors in the recoil energy 
range 5 - lOOkeV which are however consistent with expected backgrounds. For the energy 
resolution we apply (1261) as for CDMS II. We have also checked that the silicon data set 
from the same run gives no further constraints. 



3.2.5 CRESST I 

CRESST I took 1.51 kg-days of data in October 2000 [37]. The quenching factor is Q = 1. 
The data is binned from 0.6 to 20 keV, each bin covers an energy range of 0.2 keV. We 
use the events after cuts have been taken into account, cf. Figure 9 & 10 in [37]. The 
energy resolution given by the CRESST collaboration is [37] 

a(E R ) = v/(0.220 keV) 2 + (0.017£ R ) 2 . (28) 



3.2.6 CRESST II 

Two CaW04 crystal CRESST II detector modules (Verena and Zora) took data between 
March 27 and July 23, 2007. The total exposure of tungsten corresponds to 30.6 kg-days 
in an energy range of 10-100 keV with an acceptance of 0.9. Scattering off of calcium 
and oxygen can safely be neglected. Seven candidate events are shown in Figure 8 of [20]. 
For the energy resolution we us^| 

a(E R ) ~ IkeV . (29) 



There exists also a data set from an older run of CRESST II with the prototype 
detector modules Daisy and Julia [21]. Some previous studies [18,19] combined the com- 
missioning run with this older test run. We emphasize that we refrain from doing this. 
The reason is that the CRESST II detector was subject to significant modification be- 
tween the runs, in particular a neutron shield as well as a muon veto were installed [20]. 
The absence of the neutron shield in the test run in combination with the poor light 
resolution of the Julia module may have led to a "leakage" of neutron-induced recoils 
into the region of tungsten events [21]. It may well be due to this particular background 

4 The energy resolution does not exceed IkeV in the energy range up to 100 keV [45], a slight overes- 
timation of the resolution will only make our limits marginally more conservative. 
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[keVee] 




[keVee] 




0.338-0.467 


5703 


1.241-1.370 


62 


0.467-0.596 


52 


1.370-1.499 


38 


0.596-0.725 


33 


1.499-1.628 


35 


0.725-0.854 


50 


1.628-1.757 


31 


0.854-0.983 


45 


1.757-1.886 


25 


0.983-1.112 


32 


1.886-2.015 


19 


1.112-1.241 


50 







Table 3: CoGeNT binned data up to 2 keVee. 



that the observed event rate in the test run was higher than in the commissioning run. 
Therefore we find it most appropriate to include only the commissioning run in our 
analysis, which leads to stronger constraints. 



3.2.7 CoGeNT 

CoGeNT took 8.38 kg-days of data in March 2008 above energies of ~ 0.3 keVee [38, 
46]. The quenching factor is Q = 0.2. We obtained a data sample from the CoGeNT 
collaboration containing the absolute number of events in very small bins of size ~ 
30 eVee. One should, however, not use such small bins for the analysis, as gain-instabilities 
can shift the measured energy of events slightly and introduce systematical errors [46]. 
We therefore followed the procedure of the CoGeNT collaboration [38,46] and combined 
four bins into one, resulting in a bin width of ~ 130keVee, c.f. Tabled We ignored all bins 
above E' = 2 keVee corresponding to E R = 10 keV in order not to dilute the statistical 
power of the binned Poisson method. Above this energy CDMS II is very sensitive and 
CoGeNT cannot give additional limits. 

We checked that the efficiency of CoGeNT is well approximated by the formula given 
in |5], 

« (£ '»=°- 66 - 50lL- (30) 
The energy resolution used by CoGeNT is [5] 

a(E') = v/(69.7eVee) 2 + (0.98 eVee)E' . (31) 



3.2.8 TEXONO 

The TEXONO collaboration took 0.338 kg-days of data with a very-low threshold ger- 
manium detector [39] in the range E' = 0.1-7.1 keVee. The experiment was running 
between January 15 and February 6, 2007 [47]. We obtained a data sample from the 
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Figure 1: Parameter space of inelastic dark matter potentially explaining the DAMA annual 
modulation signal at 90% CL without (left panel) and with (right panel) detector resolution of 
DAMA taken into account. For the color code see text. 



TEXONO collaboration from which we extracted the number of events in bins of size 
100 eVe^E The quenching factor is Q = 0.2 and we again use only the bins up to 2keVee. 



The overall efficiency is estimated as (cf. Figure 3 in [39]) 

2 /487£' 4 

— arctan 

vr \ ke Vee 

while the energy resolution can be approximated by 



£(E') = 0.983 • 0.915 ■ - arctan ( 4 ) , (32) 



a(E') = I 0.055 + 0.009a/— — I keVee . (33) 

keVee 



4 Results 

In this section we describe the results of our analysis, i.e. we show the regions in the 
parameter space of inelastic dark matter which are compatible with the DAMA signal 
and check if they are also consistent with the null results of the various other searches 
described in Section [31 The parameter space is spanned by three parameters: the WIMP 
mass m x , the mass splitting 5 and the WIMP neutron cross section a n . 

We find three regions^ in parameter space which could potentially explain the DAMA 
signal. They are depicted in different colors in Figure CD To illustrate the impact of the 
DAMA energy resolution, the left and right panel show the DAMA allowed regions 
without and with the resolution taken into account, respectively. While the large green 



5 This data set can be obtained from the TEXONO collaboration directly. 

6 If we enhance the confidence level, the distance between the red and the green region in Figure Q] 
decreases, at the 99% level they are already slightly connected. 
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Figure 2: Recoil spectrum and modulation amplitude at DAM A for the local best fit point 
m x = 56.5 GeV, 6 = 115.4keV, a n = 2.5- 10~ 39 cm 2 . Left Panel: Differential event rate dR/dE R 
in a Nal target as a function of the recoil energy Er for the parameters above. Only iodine events 
are visible as scattering off of sodium is kinematically not accessible. The dash-dotted curve 
refers to June 2, the solid curve to December 2. Right Panel: Predicted modulation amplitude 
at DAMA for measured energies E 1 up to lOkeVee with (solid) and without (dashed) resolution 
taken into account. The modulation amplitude measured by DAMA is also included. 



and red regions are extended through the inclusion of the resolution, the small yellow 
region completely disappears. We will now discuss these three regions in turn and confront 
them with the limits from the various null searches. In this analysis we will not consider 
additional constraints e.g. from collider searches or indirect detection [48, 49] as their 
applicability strongly depends on the considered model. 



4.1 The quenched region 

The quenched iodine region corresponds to the green region in FigureEDwith vn x > 43GeV 
and < 5 < 145 keV when the detector resolution is taken into account. While the plot 
ends at m x = 1 TeV, the region extends to extremely high masses. Allowed cross sections 
range from a n ~ 6 • 10~ 42 cm 2 to 10~ 36 cm 2 , a very thin strip extends above 1CT 36 cm 2 . 

Only quenched events can contribute to the signal. This can be seen in the left panel 
of Figure [2] where the recoil spectrum in sodium iodide for one of the two local best fit 
point£| (m x = 56.5 GeV, 5 = 115.4keV and a n = 2.5- 10~ 39 cm 2 with \ 2 = 12.1) is shown. 
The lower threshold of the iodine spectra lies at recoil energies Er ~ 20 keV, i.e. there 
are no channeled events in the relevant energy range below ~ 10 keVee. 

In the right panel of Figure [2] the predicted modulation amplitude (with and without 
taking into account the detector resolution) is shown for the same parameters together 
with the measured modulation amplitude in the bins up to lOkeVee. Note that the finite 
energy resolution leads to a slight smearing of the spectrum. Especially parameter sets 
which predict too high spectra at zero resolution may become compatible with DAMA 
through this smearing. This is the reason why the DAMA allowed green region in Figured] 

7 The other local best fit point lies at S ~ 0, it has a slightly lower value of x 2 - 



16 




50 100 150 100 1000 

<5[keV] mJGeV] 





DAMA 90% 


■ 


DAMA 99% 




DAMA unmod 




CRESST II 90% 




CRESST II 95% 




XENON10 90% 




CDMS II 90% 




ZEPLIN III 90% 



Figure 3: DAMA allowed regions in the parameter space of inelastic dark matter together with 
the relevant constraints from other direct detection experiments for fixed m x = 100 GeV (left 
panel) and fixed 5 = 130 keV (right panel). The region above the lines is excluded. 



extends to larger a n and smaller m x if the resolution is taken into account. 

In Figure [3] the limits from the DAMA absolute rate and the other direct detection 
experiments on the quenched iodine region are shown. We plot only constraints from 
CRESST II, ZEPLIN III, CDMS II and XENON10, as constraints from other experi- 
ments are significantly weaker in this regime. We fix m x = 100 GeV in the left panel and 
5 = 130 keV in the right panel and scan over the two remaining parameters. For DAMA 
we show the 90% and 99% allowed regions, while all exclusion curves are given at 90% 
CL. Only for CRESST II, which gives the most important constraint, we also include the 
95% exclusion curve. In fact CRESST II excludes the whole 99% DAMA allowed region 
at 95% CL. We checked that this statement holds not only for the parameters chosen 
above, but in the complete quenched iodine region over the whole range of m x and S. 

We then tested the robustness of this statement by varying the iodine quenching fac- 
tor Qi and the local galactic escape velocity t> e sc within their experimental uncertainties. 
While high t> esc = 600 km/s and low t> esc = 500 km/s did not weaken the CRESST II 
bound, the latter becomes slightly weaker for the quenching factor at its minimal or 
maximal values of Q\ = 0.08 and Q\ = 0.1 [50] respectively. Nevertheless, even if we tune 
Qi in the most favorable way, the DAMA 99% confidence region is still excluded at 90% 
CL by CRESST II in the whole quenched iodine region. 

One should, however, be aware that we derived our limits for the Standard Halo 
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Figure 4: Same as Figure [2] but for m x = 5.82 GeV, 5 = 14.34 keV, a n = 3.9- 10" 37 cm 2 (best fit 
point in the channeled sodium region at zero resolution). The gray curves in the left panel refer 
to sodium, the black curves again to iodine. Note that the iodine spectrum has been rescaled 
by a factor 1/1000. 

Model. If the velocity distribution of WIMPs shows substantial deviations from this 
model, there may exist the possibility to reduce the level of incompatibility between 
DAMA and CRESST II. Nevertheless tungsten is kinematically favored over iodine in 
the quenched iodine regime, i.e. there are no parts in the velocity distribution of WIMPs 
which are accessible to DAMA but not to CRESST II. Consequently there is no simple 
way to resolve the tension between the two. 

4.2 The channeled regions 
4.2.1 Sodium 

The channeled sodium region which corresponds to the yellow region in the left panel 
of Figure ffl lies at 5.3 GeV < m x < 6.1 GeV, 10 keV < 5 < 16 keV and 2 • 1CT 38 cm 2 < 
cr„ < 2 ■ 10~ 36 cm 2 for the detector resolution set to zero. It is the only region where the 
DAMA modulation amplitude can potentially be explained by scatterings off of sodium. 
The reason why this region disappears when we take into account the resolution can easily 
be inferred from the left panel of Figure BJ There we show the sodium and the iodine 
recoil spectra for the best fit point at zero resolution (m x = 5.82 GeV, 5 = 14.34 keV and 
a n = 3.9-10~ 37 cm 2 with x 2 — 15.5). One can see that there is an extremely large number 
of iodine events below the DAMA threshold which in case of a nonzero resolution are 
partly shifted to higher energies. In the right panel of the same figure we see that this 
totally spoils the fit to the DAMA modulation amplitude. 

Nevertheless we hesitate to completely exclude the sodium channeled region. In some 
part of it, the iodine peak lies at energies E' < 1.5keVee. In such cases the compatibility 
with DAMA could remain if <r(1.5keVee) < 0.2keVee which is about 1/3 of the expected 
resolution of DAMA/LIBRA at that energy. One should have in mind that the resolutions 
of DAMA/Nal and DAMA/LIBRA were only extrapolated to such low energies which 
implies some uncertainties. In addition, the resolution may show some deviations from 
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Figure 5: Same as Figure OH but for m x = 5.5 GeV (left panel) and S = 13keV (right panel) and 
the detector resolution of DAMA set to zero (see text). 



Gaussianity far away from the central value. Finally there is also the possibility that the 
iodine peak is somehow affected by DAMA's noise rejection procedure near threshold. 

On the other hand, the total rate at DAMA/LIBRA (Figure 1 in [1]) shows no 
sharp rise above E' = 0.75 keVee. And while one should be very careful about the 
interpretation of bins below the threshold, there is at least no indication for the iodine 
peak down to this energy. Although the existence of the sodium channeled region is 
therefore doubtable, we nevertheless present here the limits on this region from other 
direct detection experiments. For this purpose we set the energy resolution of DAMA 
to zero. Note that we do this just for illustration, the definite position and size of the 
DAMA allowed region, if it exists, highly depends on the behavior of the DAMA detector 
below its threshold. 

As the sodium region lies at very low m x , important constraints only arise from 
the very low-threshold experiments CoGeNT, TEXONO and CRESST I. In Figure [5] we 
show the DAMA allowed regions (90% and 99%) at zero resolution together with the 
90% exclusion curves from these experiments for fixed m x = 5.5 GeV (left panel) or fixed 
5 = 13 keV (right panel). As can be seen only TEXONO has begun to explore a small 
part of the DAMA allowed region, most of it still remains unchallenged. 
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Figure 6: Same as Figure [2] but for the best fit point in the channeled iodine region (m x = 
11.9 GeV, 5 = 33.3 keV, a n = 4.7 • 10 _39 cm 2 ). All spectra refer to iodine as scattering off of 
sodium is kinematically not accessible. 

4.2.2 Iodine 

The channeled iodine region corresponds to the red region in Figure CD with 9 GeV < 
m x < 20 GeV and < 5 < 38 keV when the detector resolution is taken into account. 
Allowed cross sections range from a n ~ 9 ■ 10~ 42 cm 2 to 10~ 36 cm 2 , a very thin strip 
extends above 10 _36 cm 2 . Only channeled iodine events give a contribution to the signal. 
This can be seen in the left panel of Figure M where we plot the recoil spectrum in sodium 
iodide for the best fit point (m x = 11.9 GeV, 5 = 33.3 keV and a n = 4.7 • 10~ 39 cm 2 with 
X 2 = 12.3). As the iodine spectrum is practically zero above E R ~ 6 keV, all quenched 
events lie below the DAMA threshold and cannot contribute to the signal. 

The right panel shows the predicted modulation amplitude at the best fit point and 
the measured modulation amplitude in the bins up to lOkeVee. As in the quenched iodine 
region the inclusion of the energy resolution causes a slight smearing of the spectrum 
which leads to an extension of the DAMA allowed region. 

The most important constraints on the channeled iodine region arise from XENON10 
and CDMS II, as can be seen in Figure [7] In particular, we find that the part of the 
channeled iodine region where m x > 15 GeV is excluded by XENON10 at 99% CL. 
However, below this mass there still exists parameter space in agreement with DAMA 
and completely consistent with all other experiments. 

5 Comparison with other studies 

Our results show important deviations from earlier studies [16, 18, 19] which we now 
want to address in some detail. We focus on heavy inelastic dark matter (green region 
in Figured]) as the low mass regions have not been studied systematically. 

In contrast to [16,18,19] we took into account the energy resolution of DAMA/LIBRA 
and DAMA/Nal which leads to an extension of the DAMA allowed region especially 
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Figure 7: Same as Figure [3] but for m x = 11 GeV (left panel) and 5 = 30 keV (right panel). 



to larger cross sections and smaller masses (c.f. Figured]). We also took into account 
information from the higher energy bins. 

As for the other experiments, the most important difference concerns the constraint 
from CRESST II which is more severe in our case than in [16,18,19]. In [16] simple Poisson 
statistics were used to derive the CRESST II limit which allows a comparison between 
the absolute number of predicted and measured events only, but ignores information from 
their spectral shape. Interval based methods like the pmax-method employed in [18, 19] 
or the maximum gap method used by us seem more appropriate and lead to stronger 
limits. One would therefore expect that [18, 19] obtain comparably severe constraints. 
However, these studies combined the data set from the commissioning run of CRESST II 
with the older test run. In our analysis we refrain from including this run as explained in 
detail in Section l3.2.61 But even if this run is included there is still some tension between 
CRESST II and DAMA unless one particular event at ~ 22 keV measured by the Daisy 
module is included [19]. This event, however, enters the WIMP search region only if 
the quenching factoid of tungsten deviates from its experimental value (see discussion 
in [21]). Newer measurements of the quenching factor [51] seem to exclude this possibility, 
i.e. the mentioned event should not be treated as a candidate event. 



8 Note that in cryogenic detectors the quenching factor has no impact on the energy calibration but 
only on the discrimination between background and candidate events. 
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6 Conclusions 



In this study we examined the compatibility of the DAMA annual modulation signal 
with other direct searches for the case of inelastic dark matter. We searched for regions 
consistent with the spectral shape of the DAMA modulation amplitude by scanning over 
the three-dimensional parameter space of inelastic dark matter which is spanned by the 
WIMP mass m x , the mass splitting 5 and the WIMP neutron cross section a n . We found 
two regions in parameter space which correctly reproduce the DAMA signal, one of them 
at higher m x corresponding to quenched and one at small m x corresponding to channeled 
scattering events off of iodine. A third region where the DAMA signal can be explained 
by channeled scatterings off of sodium opens up if we decrease the energy resolution of 
DAMA. Its existence depends strongly on the behavior of the DAMA detector below 
threshold. 

All DAMA allowed regions were studied in detail and limits from the relevant other 
direct detection experiments were applied. In contrast to previous studies we found that 
for the Standard Halo Model inelastic dark matter with m x > 15GeV is ruled out at 95% 
CL where the strongest constraints are set by XENON10 (at smaller m x ) and CRESST 
II (at larger m x ). This statement remains fairly robust with respect to experimental 
and astrophysical uncertainties: when varying the galactic escape velocity and the iodine 
quenching factor within their experimental uncertainties, the DAMA 99% confidence 
region remains excluded at 90% CL. Dark matter velocity profiles which differ substan- 
tially from the Standard Halo Model may reduce the discrepancy between CRESST II 
and DAMA in the high mass region slightly. 

At lower masses, where the DAMA signal can be explained through channeled events 
at iodine and sodium respectively, we find that parts of these regions are consistent with 
all direct detection experiments. In conclusion we have shown that heavy inelastic dark 
matter is disfavored by the combination of DAMA spectral information and the exclusion 
limits from other experiments. Light inelastic dark matter on the other hand constitutes 
a viable solution to the DAMA puzzle. 
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